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Abstract 

The connection of Gribov's confinement scenario in Coulomb gauge with the 
center vortex picture of confinement is investigated. For this purpose we assume a 
vacuum wave functional which models the infrared properties of the theory and in 
particular shows strict confinement, i.e. an area law of the Wilson loop. We isolate 
the center vortex content of this wave functional by standard lattice methods and in- 
vestigate their contributions to various static propagators of the Hamilton approach 
to Yang-Mills theory in Coulomb gauge. We find that the infrared properties of 
these quantities, in particular the infrared divergence of the ghost form factor, are 
dominated by center vortices. 



1 Introduction 



The infrared sector of QCD, in particular the confinement mechanism, is not fully under- 
stood yet, although substantial progress had been made during recent years. The progress 
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comes mainly from lattice calculations [T|, which gave support to both the dual Meiss- 
ner effect [3] and the center vortex picture [3] of confinement and also indicate that these 
two pictures are likely only two sides of the same medal. Furthermore, there is lattice 
evidence that also Gribov's confinement scenario is triggered by magnetic monopoles and 
center vortex configurations [HUE]- In addition, one can show analytically that Zwanziger's 
horizon condition, i.e. an infrared diverging ghost form factor, which is at the heart of 
Gribov's confinement scenario in Coulomb gauge, implies dual Meissner effect [7]. 

In recent years there has been a renewed interest in studying Yang-Mills theory in Coulomb 
gauge, both in the continuum [8] and on the lattice [HI [101 EJ H21 H21 [13] • In particular, 
much work was devoted to a variational solution of the Yang-Mills Schrodinger equation 
in Coulomb gauge [T5l [HI [13 HH HH [201 IZH [221 [23]. Using Gaussian type of Ansatzes for 
the vacuum wave functional, a set of Dyson- Schwinger equations for the gluon and ghost 
propagators was derived by minimizing the vacuum energy density. An infrared analysis 
|22| of these equations exhibits solutions in accord with Gribov(-Zwanziger) confinement 
scenario. Imposing Zwanziger's horizon condition one finds an infrared diverging gluon 
energy and a linear rising static quark (Coulomb) potential [22] and also a perimeter law 
for the 't Hooft loop |23| . These infrared properties are reproduced by a full numerical 
solution of the Dyson- Schwinger equations over the entire momentum regime [IH] , and are 
also supported by lattice calculations jlO ^ [T3 ] [T^] . Moreover, these lattice calculations also 
show clear evidence for an infrared divergent ghost form factor and an infrared suppressed 
static gluon propagatoiQ. While the horizon condition has to be imposed by hand in 
D — 3 + 1, as there exist also subcritical solutions with an infrared finite ghost form 
factocl |2U) . the coupled set of Dyson-Schwinger equations in D = 2 + 1 only allows for 
critical solutions having an infrared diverging ghost form factor [23]. Furthermore, recent 
studies within the functional renormalization group treatment of the Hamilton approach 
in Coulomb gauge yields the horizon condition as a solution of the flow equation [25] . In 
D = 1 + 1, finally, the exact ghost form factor is infrared enhanced 



lr rhis is different from Landau gauge, where lattice calculations seem to indicate an infrared finite 
gluon propagator and ghost dressing function [28], in contradiction to the scaling solution of the Dyson- 
Schwinger equations |35) . 

2 These subcritical solutions are the analogue of the so-called de-coupling solution of the Dyson- 
Schwinger equations in Landau gauge. 

3 Yang-Mills theory in D = 1 + 1 is non-trivial only on a compact manifold. In the Hamiltonian 

approach (with a continuous time) space is Si so that the momenta are discrete. The zero momentum 
corresponds to a zero mode of the Faddeev- Popov kernel and is therefore excluded [21 j . 
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Despite the encouraging results for the infrared properties of the various Greens function 
(in particular the linear Coulomb potential) and the good agreement of the gluon energy 
|2U] with the lattice data [131 E], the crucial test for the wave functional in the variational 
approach, namely the calculation of the Wilson loop, still has to come. In fact, a linear 
Coulomb potential is necessary but not sufficient for confinement since the Coulomb string 
tension is only an upper bound to the Wilsonian string tension [S]. With the variational 
vacuum wave functional at hand it would be, in principle, straightforward to calculate 
the Wilson loop. However, path ordering makes an exact evaluation of the Wilson loop 
impossible. In Ref. [26J the spatial Wilson loop was calculated from a Dyson equation, 
which takes care of the path ordering in an approximate fashion, employing the static 
gluon propagator as input. Despite the rather limited range of applicability of the Dyson 
equation, a linearly rising potential could be extracted from the obtained Wilson loop. 
Here we will proceed along a different line. We assume a wave functional which is known 
a priori to produce an area law for the (spatial) Wilson loop, and calculate with this wave 
functional the various propagators of the Hamiltonian approach in Coulomb gauge. The 
infrared properties of these propagators are then compared with the ones obtained from 
lattice results in D = 4 Coulomb gauge which, in turn, agree qualitatively with findings 
from a Gaussian type of variational wave functional. 

A simple choice for a confining wave functional is 



where denotes the spatial components of the non-Abelian field strength and \i is a 
dimensionfull parameter, which, in principle, could serve as variational parameter. We 
will later see that the scaling properties of D = 3 Yang-Mills theory tie the parameter \i 
to the numerical value of the (spatial) string tension: \x thus merely sets the overall scale. 

The wave functional in eq. ([1]) models the infrared sector of the Yang-Mills vacuum: it 
is gauge invariant and can be considered as the leading order in a gradient expansion of 
the true Yang-Mills vacuum wave functional (27]. Furthermore, the functional of eq. (00) 
produces an area law for the (spatial) Wilson loop. This is because the D = 4 expectation 
value of any gauge invariant and Ao-independent observable Q[A] in this state is precisely 
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given by the one in the D = 3 dimensional Yang-Mills theory, 




(2) 



For gauge variant observables such as the Green functions, we need to pick a specific 
gauge on the vector potential A. Choosing Coulomb gauge VA = 0inD = 4 obviously 
entails Landau gauge for the D = 3 YM-theory. 

The wave functional of eq. ([Q) is certainly inappropriate at large momenta where it yields 
a gluon energy cu(|k|) ~ |k| 2 instead of u(k) ~ |k|. However, this should be irrelevant 
for the confining properties. Also in the deep IR region one does not expect eq. (CQ) to 
exactly reproduce the D = 4 Yang-Mills Theory, since the standard lattice gluon and 
ghost propagator in D = 3 Landau gauge rather satisfy a decoupling type of solution 
|28| . in contrast to D = 4 Coulomb gauge [T3| IT4"] . Indeed, one would expect the correct 
wave-function to be better described by |29) : 




where D 2 is the adjoint covariant Laplace operator. In Ref. |29| . where the D = 2 + 1 
theory was examined, the choice c = — Ao + m 2 was made, Ao being the lowest eigenvalue 
of — D 2 . The technical difficulties of working with eq. ([3]) are however beyond the scope 
of this paper, since one does not expect the Laplacian term to modify the vortex content 
of the D = 3 theory. Even the simplified version given in eq. (OQ) cannot be used for 
analytic calculations; instead, the calculation of expectation values in this state requires 
conventional 3-dimensional Yang-Mills lattice simulations. As stated earlier, these D = 3 
lattice calculations must be carried out in Landau gauge when ([T]) is used as Schrodinger 
wave functional in Coulomb gauge. 

The infrared properties of the vacuum sector of Yang-Mills theory are known to be dom- 
inated by center vortices [I]. This is also true in D = 3. The wave functional eq. ([T|) thus 
contains, among other configurations, an ensemble of percolating center vortices. By stan- 
dard lattice methods [30] we can extract the center vortex content of this wave functional. 
Therefore the use of the wave functional eq. ([I]) also allows us to study (on the lattice) 
how, within the Hamiltonian approach, the Gribov-Zwanziger confinement mechanism is 
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related to the center vortex picture of confinement. Previous lattice calculations jH |6] 
have shown that removal of center vortices from the (4-dimensional) Yang-Mills ensemble 
by the method of Ref. [30] makes the ghost form factor infrared finite, analogously to 
the suppression observed before in Landau gauge pD] . In the present work we will calcu- 
late the center vortex contribution to various Coulomb gauge propagators. The paper is 
organized as follows: 

In the next section, we briefly recall some properties of D = 3 Yang-Mills theory, in par- 
ticular the scaling behavior and some known facts about Landau gauge Green functions. 
Section 3 contains a description of our numerical setup and the results for the static gluon 
and ghost form factors in Coulomb gauge. We also discuss the role of center vortices and 
their implications for the Gribov-Zwanziger scenario. We close with a brief summary and 
an outlook on future investigations. 



2 Yang-Mills theory in three dimensions 

Since we want to describe the continuum model of eq. ([I]) using a lattice, we must first 
have a closer look at the scaling properties of D = 3 YM theory. The lattice model is 
defined on a D = 3 cubic space-time grid with periodic boundary conditions. We will 
use the Wilson action and employ various gauge fixing algorithms. The scale fi of the 
continuum wave functional in eq. ([T|) plays the role of the three-dimensional bare YM 
coupling in the continuum limit, \i = g\. 

In three-dimensional YM theory, the renormalized coupling g^ is 

a 9o 

in terms of the bare coupling g and the lattice spacing a. The string tension in units of 
the lattice spacing is therefore 

9 16 16 cr 

Since gn sets the overall scale, cr/g^ is a dimensionless constant independent of (3, i.e. large 
Creutz ratios should scale according to 

X (R, R) « a((3) = ^ + ^(/T 3 ) • (5) 
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Large-scale simulations in SU(2) |32| have indeed confirmed the existence of a scaling 
window (3 ~ 3 ... f 2 in which the dependency can be observed. More precisely, the 
best fit |32] is 

1.788 ( . 1.414 \ 
valid for /3 > 3. From eqs. (j4]), ([6]) and \i = g$, we find 

' . (7) 



a 

— = 0.111 

ft 2 



1.414 

1 + — + 



i.e. the variation parameter \x = fi(a, (3) is fixed once the overall scale a and the lattice 
coupling (3 in the scaling window are given. For the standard value <Jo = 440 MeV, fi 
falls in the range \i = 1088 MeV . . . 1236 MeV when /3 is varied in the scaling window. 
There is however no compelling reason for the string tension of this D = 3 model of the 
D = 4 Yang-Mills vacuum to coincide with the "physical" string tension of the genuine 
D — 3 Yang-Mills theory. In view of the following we will then treat \x as an adjustable 
parameter that controls the tradeoff between matching the D = 4 string tension or the 
corresponding Green functions. 

Ghost and gluon Green functions in Landau gauge are generally expected to be multi- 
plicatively renormalizable, i.e. the lattice data for different /3 must fall on top of each 
other once the momentum is expressed in physical units and a finite field renormalization 
Z(/3) is applied. The Gribov-Zwanziger confinement criterion makes qualitative predic- 
tions about the deep IR behavior of Green functions. It is based on the idea that the 
gauge field configurations in the functional integral should be restricted to the so-called 
fundamental modular region (FMR) . The dominant contributions within the FMR should 
then come from field configurations near the Gribov horizon, where the near-zero modes 
of the Faddeev-Popov operatoiQ M = —3D strongly enhance the Coulomb potential 



K(|x-y|)~Tr (M-^-AJM- 1 ). 



To enforce the restriction to the FMR, a non-local but renormalizable horizon term may 
be added to the Yang-Mills action. The fact that the partition function is dominated by 
near-horizon configurations is then expressed as the so-called horizon condition, which in 
turn implies that the ghost propagator should be more singular than the free ghost in the 



We write D = 9 M + g A — — gf ' A c ^ for the covariant derivative acting on adjoint color fields. 
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deep infrared, 



G a \p) 



sab 

((M~ 1 ) ab ) — ~-rd(p) 



lim d~\p) = 0. 



(9) 



Similarly, the gluon propagator should be infrared suppressed or vanishing, 



(p)> = ** 8, 




) 



D{p) 



lim D[p) < oo. 



(10) 



These results are expected to hold both in D = 3 and D = 4 dimensions. It should be 
noted, however, that a soft BRS breaking but renormalizable mass term is possible in 
D = 3 and D = 4 which gives rise to the so-called decoupling solution. The latter is 
characterized by a tree-level-like ghost propagator in the IR (<i(0) = const) accompanied 
by an infrared finite, but non-vanishing gluon propagator -D(O) ^ 0. Recent lattice 
investigations seem to favor this type of behavior in Landau gauge for both D = 3 and 
D = 4 which is, however, at odds with both the Gribov-Zwanziger and the Kugo-Ojima 
criterion favored by functional methods in the continuum. As for the Coulomb gauge 
in D = 4, there is, however, ample evidence both from variational methods fT7\ IT9] and 
lattice investigations [T3| IT4"] that the ghost propagator is infrared enhanced while the 
static gluon propagator is infrared vanishing, in accord with the original Gribov scenario. 
In the next section, we compute the gluon and ghost Green functions in the deep infrared, 
using the confining wave functional eq. flTJ as a model of the D = 4 ground state in 
Coulomb gauge, reliable at least in some intermediate momentum range. Technically, 
this amounts to a lattice simulation in D = 3 Landau gauge; in contrast to earlier studies 
cited above, we will first perform a detour via the maximal center gauge (MCG) to identify 
the percolating vortex content of each configuration. We can then remove or isolate these 
vortices before going to Landau gauge, and thus study the interplay between the vortex 
and Gribov-Zwanziger scenario. 

3 Numerical results 
3.1 Numerical setup 

Since we are mainly interested in the IR properties we choose in the following a fixed 
coupling (3 = 3.5 close to the smaller end of the scaling window, which in turn gives access 
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to the smallest momenta while still allowing to extract continuum physics. Calculations 
were carried out on 20 3 and 40 3 lattices using a standard over-relaxed version of Creutz' 
heath-bath algorithm. We employed 500 sweeps for thermalization and 50 to 100 sweeps 
between measurements to reduce auto-correlations. For each Green function, we took 300 
measurements on the smaller volume 20 3 , and 100 measurements for the volume 40 3 . 
In each measurement, the lattice configuration was first brought to maximal center gauge 
(MCG) to allow for a clean identification of the center vortex content. The vortices were 
then either removed or isolated, and both fields (in addition to the unmodified MCG 
configuration) were subjected to further Landau gauge fixinglf] As explained above, this 
corresponds to D = 4 simulations with the trial wave functional eq. (including its 
center vortex content) in Coulomb gauge. 

For both gauge fixing steps, we used an iterated over-relaxation algorithm with up to 30 
restarts after a random gauge transformation on the initial configuration to reduce the 
Gribov noise. The gauge fixing procedure was terminated when the local gauge violation 
was sufficiently small, 

HAIL = max max \6 k (x) I < e = 10~ 12 , 

where in terms of the the quaternion representation 11^ = a? + ia k of the SU (2) links, 

1 2 

MCG : 6\x) = - Wx) a%x) - a^x - jl) a%x - /*)] 

At=0 
1 2 

Landau : d\x) = - ^ (aj(x) - aj(x - fi)) . 

In all cases, the gauge fixing functional was numerically stationary long before the viola- 
tion limit was reached. 

For Landau gauge, we further improved the quality of the gauge fixing using either flip 

pre-conditioning steps |33] and, on the larger 40 3 lattice, simulated annealing methods. 

However, these improvements had only a marginal effect on the attained maximum of the 

gauge fixing functional, since the previous MCG fixing acts as a strong pre-conditioning 

even when vortices were not removed. 

5 It should be noted that MCG corresponds to Landau gauge for the adjoint link, which reduces to 
ordinary Landau gauge once the singular vortex background is removed. Thus, the vortex removed MCG 
configurations are already in Landau gauge. Nevertheless, we have carried out the subsequent Landau 
gauge fixing step to be sure that we converged to the best minimum within our numerical precision. 
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Figure 1: Left panel: The D = 3 gluon propagator in Landau gauge after prior MCG fixing, 
both with and without vortices. Right panel: The same plot for the D = 3 ghost form factor in 
Landau gauge. 

3.2 Gluon propagator 

The gauge field on the lattice is extracted using the standard definition, 

a((3) A%x + A/2) « i [U,{x) - Ul(xj) = 2a» 

where a((3) is the lattice spacing and = a° + ia is the quaternion representation 
of the links§ A fast Fourier transformation then gives access to the gluon propagator 
eq. (HQJ0 

In the left panel of Fig. [TJ we first show the D = 3 gluon propagator in Landau gauge, 
after prior MCG fixing. As long as vortices are not removed, the resulting propagator is 
virtually identical with the one obtained by fixing minimal Landau gauge directly. As can 
be clearly seen, the removal of vortices has little effect in the UV but clearly suppresses 
the propagator at intermediate and small momenta. 

Fig. [2] shows the equal-time gluon propagator Dc(p) for Coulomb gauge in D = 4, as 
obtained before removal of scaling violations. We have compared this data to the D = 3 
gluon propagator 

6 This assumes that the gauge fields are sufficiently smooth and the links close to unity, a° « 1, as 
ensured by Landau gauge and the continuum limit. 

7 Since only a single coupling /3 = 3.5 was considered, no further renormalization was necessary. On the 
smaller lattice, we have verified that all considered ensembles do lead to multiplicatively renormalizable 
gluon and ghost propagators. 
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Figure 2: The equal-time gluon propagator in D = 4 Coulomb gauge, compared to the D = 3 
gluon propagator in MCG+Landau gauge, both with and without vortices. In the left panel, 
we have only included D = 4 data from two close couplings, which hides the inherent scaling 
violations. If more couplings are included (right panel), the scaling problems become apparent 
and the comparison with the D = 3 data is less favorable. 

a. directly in (minimal) Landau gauge 

b. in Landau gauge after prior MCG fixing, 

c. in Landau gauge after prior MCG fixing and vortex removal 

(The cases a. and b. are virtually identical, and only b. is shown in Fig. [2] to make 
the plot less cluttered.) The confining D = 3 MCG/Landau propagators -Dl(p) a. and 
b. show good agreement in the low and intermediate momentum range with the naive 
instantaneous D = 4 Coulomb propagator, once the (Wilson) string tensions in the two 
calculations are set equal@ (We took yfa = 440 MeV to set the scale on the momentum 
axis.) As expected, the non-confining Landau gauge propagator with vortices removed 
cannot be matched with the D = 4 Coulomb data in any momentum range. 
However, this agreement of the confining propagators must be considered spurious, since 
the naive D = 4 Coulomb data is subject to severe scaling violations. In Ref. |13| . these 
violations were thoroughly analyzed and a procedure to extract the true instantaneous 
propagator Dq{p) in the infinite volume and continuum limit was devised. The result 

8 In the UV, we have the expected deviations, since the 1/p behavior of the Coulomb result decays 
much slower than the 1/p 2 typical for the Landau case. 
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Figure 3: The equal-time gluon propagator in D = 4 Coulomb gauge, compared to the D = 3 
gluon propagator in MCG+Landau gauge, both with and without vortices. The scaling violations 
in the Coulomb data are removed |13| to exhibit the true continuum result. In the left panel, the 
(Wilson) string tension of the D = 3 and D = 4 calculations are matched, while the right panel 
has 0"3£) « 6.2o"4£). (For simplicity, we have scaled cj^d and left a^D at the standard value.) 

displays proper scaling and the correct deep ultra-violet behavior Dc(p) ~ IpI" 1 ; it can 
also be fitted in the entire momentum range through the Gribov formula |13j . 
As can be seen from the left panel in Fig. [31 the improved static Coulomb propagator 
-Dc(p) differs considerably from the naive result, so that the spurious agreement with the 
confining D = 3 ensembles is destroyed. In particular, the deviation between the strong 
1/p 2 decay of the D = 3 propagator and the slow l/|p| decay of Dc(p) is now much more 
apparent. Qualitative agreement in the phenomeno logically important momentum range 
around 1 GeV can only be reached if the D = 3 momentum is scaled with a factor ~ 2.5 
as compared to the D = 4 case (see right panel in fig. [3]). This in turn would mean that 
the scale parameter fj, in the initial wave functional eq. ([I]) must be altered by the same 
factor and the D = 3 string tension would no longer match the D = 4 result Jj The best 
agreement is thus obtained if a^o ~ 6.2 ct^d, corresponding to \i ~ 7 GeV. 
The vortex-removed ensemble d. leads to a gluon propagator which is incompatible 
with the D = 4 Coulomb result in any momentum range, even if the restrictions on the 
parameter /i are relaxed. This shows that the percolating center vortex content is not only 
indispensable for the confining properties of the Wilson loop, but also for the qualitative 



9 In the very deep IR, the D = 3 and D = 4 propagators cannot be matched even with the relaxed 
condition on /i, since the correct D = 4 static Coulomb propagator vanishes as |p| — > 0, while the D = 3 
MCG-Landau propagator D^(p) goes to a non-zero value in this limit. 
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x D=4 Coulomb (scaling improved) 

o D=3 MCG + Landau 

+ D=3 MCG + Landau [vortex removed] 
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Figure 4: The gluon static propagator in D = 4 Coulomb gauge, compared to the approximation 
pDi(p) of eq. (|lip in D = 3 MCG+Landau gauge. 

behavior of the gluon propagator at intermediate momenta. 

As stated above, a better approximation of the correct vacuum wave functional is given 
by eq. ([3]). Indeed, either by inspecting Fig. [3] or on simple dimensional arguments, it 
should be clear that one should rather compare the D = 4 static Coulomb propagator 
Dc{p) with pDl(p) [13], which on the other hand will mimic the somewhat less crude 
approximation of the vacuum wave functional: 

tf[A] = A/-ex P (~ J tfxFvW-^FvW) , (11) 

The data basically coincide when tuning er 3 £) « 1.49 cr^, as can be seen in fig. HI Such 
good agreement in the direction of the wave functional given in eq. Q) definitely deserves 
further investigation. 
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Figure 5: The ghost form factor in D = 4 Coulomb gauge, compared to the same quantity in 
D = 3 MCG+Landau gauge. The D = 4 data is from a 32 4 lattice at j3 = 2.2 and /3 = 2.3, while 
the D = 3 data was obtained on 40 3 lattices with f3 = 3.5. 

3.3 The ghost form factor 

Although large volume simulations performed in a setup similar to ours clearly show a 
decoupling behavior also for the Landau ghost in D = 3 [28], for the momenta avail- 
able in our simulations this regime still has not kicked in, so that eq. (CQ) might still 
be considered a valid approximation in the following. Future simulations closer to the 
thermodynamic/continuum limit will however explicitely have to deal with eq. ([3]) or at 
least eq. (|TT|) . since the simple dimensional argument leading for the static gluon to the 
comparison in fig. H] cannot obviously be applied to the Faddeev-Popov operator. 
The right panel of figure [T] shows the dimensionless ghost form factor (jH]) as obtained 
from D — 3 Landau gauge with prior MCG fixing. Without vortex removal, the result 
is virtually identical to the direct Landau gauge fixing, and a power law behavior of the 
form 

d(p)~ p - 2K K = 0.23(1) (12) 
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can be fitted. This agrees well with the infrared exponent of D = 4 Coulomb gauge 
|10) . which is not subject to scaling violations. Fig. compares the ghost form factors 
in D = 3 Landau gauge and D = 4 Coulomb gauge quantitatively; as expected from the 
matching IR exponent, the comparison is fairly good. In the UV, the d(p) approaches 1 
in both D = 3 and D = 4, so the qualitative agreement is also good in this regime. (A 
closer comparison of possible anomalous dimensions in the sub-leading terms is beyond 
our current numerical precision.) 

It should be noted that the power-like behavior eq. f|T2|) is conformally invariant, i.e. the 
exponent k is not affected by any rescaling of the momenta or string tension. This allows 
to maintain the good agreement of the ghost form factor, while relaxing the condition on 
fi to make the gluon propagators match (approximately) . With the higher value of \x and 
the 3D string tension, it is therefore possible to reproduce at least qualitatively both the 
gluon and ghost form factor from the wave functional eq. ([T|), at least in the intermediate 
momentum range. 

From the right panel in Fig. (TJ it is also seen that the vortex-removed ghost form factor 
in D — 3 is virtually flat (k < 0.05) so that it cannot be matched even qualitatively 
to the D = 4 Coulomb result. In particular, it seems impossible to satisfy the horizon 
condition d(O)" 1 = once center vortices are removed. This is, of course, expected since 
vortices live on the Gribov horizon and are thus indispensable to maintain enough near- 
zero modes of the Faddeev-Popov operator to produce the diverging ghost form factor. 
A similar conclusion was also reached in Ref. [6] from the direct study of the low-lying 
Faddeev-Popov spectrum. 

4 Conclusion and outlook 

In this paper we have investigated the properties of the static propagators of the Hamilto- 
nian approach in Coulomb gauge, assuming the wave functional eq. ([ID, which shows strict 
confinement, as an approximation of the low energy limit of the true Yang-Mills vacuum 
wave functional. The calculation of expectation values of observables in the D = 3 + 1 
Hamiltonian approach in Coulomb gauge in this state requires and ordinary D = 3 lattice 
simulation in Landau gauge. Adjusting the free scale fi of the state eq. ([I]) yields propa- 
gators which in a low to intermediate momentum regime (up to about 1 GeV) reproduce 
quite well the exact lattice propagators of the D = 3 + 1 Coulomb gauge, which also 
agree qualitatively with propagators obtained in the variational approach to continuum 
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Yang-Mills theory in Coulomb gauge. An even better agreement should be obtained by 
employing the wave functional given in eq. (J3|), (|lip . as Fig. H] shows. 
The propagators drastically change in the infra-red when the magnetic center vortices 
occurring with the weig ht |*[A]| 2 in the wave functional are removed by the method of 
Ref. [30J . In particular, the ghost form factor loses its infrared singularity and the horizon 
condition is no longer satisfied. These results indicate that there is no need to include 
additional explicit center vortex degrees of freedom in the trial wave wave functional of 
the variational approach, as was proposed recently |34] . 
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